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Gravitational and mass distribution effects on stationary 
superwinds II. Extended dark matter haloes 
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ABSTRACT 

In this second part, we generalize the results of the previous paper We present an analytic 
superwind solution considering extended gravitationally-interacting dark-matter and baryonic 
haloes. The incorporation of the latter is critical, since they can have a substantial effect on the 
hydrodynamics of superwinds generated by massive galaxies. Although the presence of ex- 
tended and massive haloes does not change the limit for the closed-box enrichment of galaxies 
established in the first paper, they can trigger an earlier activation of the open-box enrichment 
scenario, since their gravitational potentials can contribute to the inhibition of the free su- 
perwind. Moreover, the incorporation of the extended haloes will also enhance the physical 
setting behind the superwind model, as we consider mass distributions with properties that 
emulate the results of recent simulations of ACDM haloes. 
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1 INTRODUCTION 

In the first part of this series of papers (Anorve-Zeferino & 
Corona-Galindo 2010; hereafter Paper I), we presented a simple, 
spherically-symmetric galactic superwind model considering non- 
uniform dynamical mass distributions with analogue energy and 
mass injection rates. Galaxies were modeled in terms of four pa- 
rameters: a characteristic object radius, rsc; the effective energy 
deposition rate, Ecs\ the effective mass deposition rate, i\fcff ; and 
a normalized spatial distribution, ps- The latter defined the distribu- 
tions of the dynamical mass and of the mass and energy injection 
rates within r < rac- 

The spatial distribution was assumed to follow a truncated ver- 
sion of the Dehnen profile (1993), which allows to recover trun- 
cated versions of a plateau-like and the Hemquist (1990) and Jaffe 
(1983) profiles as particular cases. As an initial simplification, we 
considered only the dynamical mass contained within rsc, which 
was assumed to account for most of the galaxy dark matter (DM) 
and baryonic mass (BM). This simplification allowed us to make 
a direct comparison between our analytic formulation and the nu- 
merical results of Silich et al. (2010), whom assumed a uniform 
distribution of the relevant galaxy parameters. Since in their model 
they also considered only the gravitational field of the central re- 
gion, we were able to analytically reproduce their numerical results 
as a particular case. However, the previous assumption is a maxi- 
mal extrapolation of the results of Persic, Salucci & Stel (1996) and 
Salucci & Persic (1997), whom found that in some cases the pres- 
ence of dark matter begins to be important well within the galaxy 
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optical radius, Rop%- They estimated that from the galaxy centres 
up to i?opt, the fraction of DM goes from 0% up to 30%-70% (see 
also Salucci et al. 2007). 

Thus, the initial supposition, as well as the assumption of 
spherical symmetry, can only be adequately interpreted in the con- 
text of a zeroth-order approximation for evaluating the effect of 
the gravitational field of ellipsoidal galaxies on the inner super- 
wind solution (r ^ rsc)- Similarly, the zeroth-order distortion of 
the external superwind hydrodynamical profiles (r > rsc) is only 
adequately predicted for galaxies with either low masses or very 
diffuse DM and BM haloes. 

For disc galaxies with large masses, the aforementioned sim- 
plification certainly does not hold, as the observed flat rotation 
curves of the extended discs require a significant amount of dark 
and baryonic matter outside of the bulge (rsc ~ '"bulge)- Needless 
to say, the presence of discs will also produce collimated flows. 
Furthermore, extended and fairly massive haloes can also have im- 
portant repercussions on the superwind solution for the case of el- 
lipsoidal galaxies (rsc ^ J'nucicus), according to the properties of 
DM haloes derived by Persic et al. (1996) and Salucci & Persic 
(1997). 

Nevertheless, the zeroth-order superwind solution that we pre- 
sented in Paper I has the advantage of being analytic. Thus, it can be 
used to construct the solution for the case in which departures form 
spherical symmetry are important. In order to do this adequately 
for a wide range of galaxy masses, we need first to incorporate the 
effect of the external DM and BM halo under the assumption of 
spherical symmetry. 

In this work, we obtain such a solution considering 
gravitationally-interacting external haloes (Section[2j. We consider 
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halo distributions that emulate the results of recent cosmological 
simulations (Section[3]l- Later, in Section|4l we obtain new limits 
for the open-box enrichment scenario and the existence of acceler- 
ating superwind solutions (see Paper I). Finally, in Section [5] we 
evaluate the impact of the halo gravitational potential on the su- 
perwind hydrodynamical profiles. The conclusions are presented 
in Section[6] 



2 AN ANALYTIC SOLUTION FOR THE FREE 
SUPERWIND INCLUDING EXTENDED HALOES 

When the haloes are included, the equation of conservation of en- 
ergy outside of the galaxy characteristic radius, equation (9) in Pa- 
per I, transforms into 



1 d[pKr^(ln^ + (77 + l)£ 

J.2 jj^ 



-/9m(V<^ + V(?!)h) 



(1) 



Above, the hydrodynamical variables are represented by 
their usual symbols, rj is the polytropic index and — Vt^i = 
—GMom/r , where Mom is the total dynamical mass within rsc. 
Similarly, — V(j!)h = —GM-ti{r)/r^, where Afh(r) is the cumula- 
tive dynamical mass (i.e. DM-l-BM) of the external halo, which has 
a total mass Mh. We will allow the profile of the external halo to 
be defined either as a continuation of the internal profile or as a 
centrally truncated profile with different characteristics. 

The integration of equation ^ yields a Bernoulli-like equa- 
tion 



lu^ + (ri+l)^ = ^V^ -TCP-, 



(2) 



where </> — —GMoM/r, (j}h{r) is the gravitational potential at 
r > rsc associated to the non-truncated version of the external 
halo, T = 1 — Mh /Mdm accounts for truncation effects, and Vg is 
the asymptotic terminal speed, which is given by 
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(3) 



In the last equation, Voo is the effective terminal speed due to 
the thermalization of SNe ejecta and individual stellar winds inside 
of the galaxy, a determines the steepness of the central dynami- 
cal mass distribution (r < r^c) and A its concentration [see equa- 
tion (11) in Paper I]. Note that r = implies an uninterrupted, 
continuous gravitational potential. Similarly, r < implies a cen- 
trally truncated external halo with a mass larger than Mum, and 
< r < 1 implies the opposite. When r = 1 there is no external 
halo, and thus (ph is identically zero. 

As in Paper I, we will work in terms of dimensionless vari- 
ables. For the present case they are: 



where T4g is given by 

_ _ 2GMriM 

V cr 
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and $h(i?) is (j>h{r) written in terms of R and normalized to Vg /2. 
The conservation laws can now be reduced to the same governing 
differential equation than in Paper I, see its equation (43). 

Thus, within the theoretical framework developed in Paper I, 
it is very easy to prove that the transonic free superwind solution is 
given by 
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Again, as in our previous work, we will give preference to the 
parametric version of the solution: 

x[x-TV,g+Rii^{R)f^ = Do{y+iy^ {2iq-y)-^ ,{■&) 

where j/ is a parameter that varies between and 2r; and 

a: = i? + rKg--R$h(i?), (9) 



y = 
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(11) 



To obtain the hydrodynamical profiles, one just needs to fol- 
low the algorithm presented at the end of section 3.3 in Paper I. An 
advantage of the parametric solution is that it allows to work with 
just functions of R in the first two critical steps, related to equa- 
tions {S]) and (|9}. On the other hand, equation ^ involves both R 
and U . For r; = 3/2 (equivalent to the case 7 = 5/3, where 7 is the 
adiabatic index) there is no need for a numerical root finder in the 
first step of our algorithm. In the second step however, its use will 
be most likely unavoidable, as the particular form of the assumed 
gravitational potential (i.e. of the external halo profile) is involved. 

In Section|4l we will give the limit above which the stationary 
solution is disrupted in the external zone (r > rsc) and the neces- 
sary condition for an accelerating stationary superwind solution. In 
order to do this, we will specify first the normalized potential c&h 
in the next section. 



3 THE EXTENDED HALO PROFILES 

How are the DM and BM distributee^ outside of the galaxy charac- 
teristic radius? Since we have permitted centrally truncated profiles 
for the external halo, theoretically, we can choose practically any 
of the usually assumed distributions; e.g. a NFW profile, Navarro, 
Frenk & White 1997; a generalized NFW profile, Moore et al. 
1999; a self-similar profile, Yoshikawa & Suto 1999; an isother- 
mal profile, and so on. Given that the most commonly used profiles 
depend on at least two parameters, and given also the additional 
freedom introduced by our truncated halo scheme; there is a vast 
number of profiles and parameters that can give reasonable agree- 
ment with observational studies and with the predictions of cosmo- 
logical simulations. 

We will try to rely on physical insight for selecting the exter- 
nal halo profile that we will use in our model. Recent cosmological 
simulations carried out by Abadi et al. (2010) predict that dark mat- 
ter haloes always contract as a result of galaxy formation. They also 
found that the contraction effect is substantially less pronounced 
than predicted by the adiabatic contraction model (Blumenthal et 
al. 1986). On similar grounds, according to the high-resolution N- 
body cosmological simulations of ACDM haloes carried out by 

^ As in Paper I, we will assume that their distributions are analogue. 
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Navarro et al. (2010), the departures from similarity in the veloc- 
ity dispersion and density profiles correlate in such a way, that a 
power law for the spherically averaged pseudo-phase-space den- 



sity is preserved, p/cr oc r 



They remarked that the index 



of the previous power law is identical to that of a Bertschinger's 
similarity solution for self-similar infall onto a point mass (in an 
Einstein-de Sitter Universe). They conclude that ACDM haloes are 
not strictly universal, but that the departure from similarity previ- 
ously mentioned may be a fundamental structural property. 

Bearing in mind the results described above, we deduce that 
the cases r < and < r < 1 correspond to artificial mathe- 
matically induced constraints that make continuous the potential at 
r = Vac for arbitrarily-chosen external-halo profiles [see equation 

The case r = 1 is physical, but corresponds to the zeroth- 
order approximation regarded as inadequate for some galaxies in 
Section[T] Thus, the case r = is of special interest, as it implies 
an unforced continuity of the gravitational potential. It turns out 
that adequately chosen truncated Dehnen profiles satisfy naturally 
the latter condition. 

For r < Tsc, the cumulative dynamical mass corresponding to 
a truncated Dehnen profiltQ is given by equation (13) in Paper I: 



M{r)^MoM{l + Af 



R 



R + A 



(12) 



We will also assume a truncated Dehnen profile for the exter- 
nal halo, but we will demand a cumulative mass of the form: 



Mh(r) = Mdm(1 + Al) 



R 



R + Ai 



3 — ai 



(13) 



At i? = 1 we have that M(l) = Mh(l). Note that for this, 
we do not require Ai = A nor qi — a. The last property can 
be interpreted in terms of a contraction of an initial spatial con- 
figuration of DM and BM with concentration Ai and steepness 
ai which produced a new configuration with concentration A and 
steepness a for r < Vsc, or well, vice-versa, if other processes 
were involved (v.gr. angular momentum). On the other hand, a triv- 
ial but important relationship can be obtained from the condition 
M(l) = Mh(l) by separating the baryonic and dark matter com- 
ponents: 



(14) 



This could be interpreted as an integral equivalent of the equation 
for adiabatic collapse derived by Blumenthal et al. (1986). Addi- 
tionally, given that the radial velocity dispersion associated to the 
Dehnen profile goes as cr ~ r"^^ when r — >■ 0, we are able to 
recover the index of the Bertschinger's power law near the centre 
of the galaxy when a = 3/4. However, Navarro et al. (2010) ob- 
tained the index from radial averaging, which implies that a can 
adopt values within a wider range. 

Note that in turn, the previous configurations could be inter- 
preted as the result of the contraction of an unperturbed configura- 
tion away from the galaxy. This is equivalent to saying that a galaxy 
formed from the perturbation of an initial state (AcQo), and that 
after certain time, the perturbation bifurcated and produced two 
inner contracted states characterized by (A,a) and (Ai,ai). The 
first state characterizes the inner regions of the galaxy, r < Vsc- 
Then, the characteristic radius Vsc can be taken either as the ra- 
dius of a galaxy nucleus or of a bulge. The second state charac- 
terizes the outer portions of the galaxy (e.g. a disc -I- DM). This 



is in agreement with the aforementioned cosmological simulations, 
and it implies that galaxies carved out gravitational potential holes 
when they formed, and that they correspond to local depressions of 
an otherwise smoother gravitational potential. 

Here, we are just interested in the superwind solution, so, in 
order to keep things simple, we will just consider the states {A, a) 
and (Al, ai), i.e. we will ignore the depression of the reference 
gravitational potential (y4o,ao). The price that we will pay for 
this, as well as for the implicit analogue distribution of the bary- 
onic and DM components assumed in our scheme, is that instead 
of (almost) 'perfectly' flat rotation curves up to 15 times the opti- 
cal radius (Persic et al. 1996, Salucci & Persic 1997), the rotation 
curves will show some downwards skewness at large radii. They 
are however well above the curves corresponding to keplerian ro- 
tation of the baryonic mass. Evenmore, the behaviour of the asso- 
ciated rotation curves away from Vsc is consistent with that of the 
universal rotation curves derived by Salucci et al. (2007) for spiral 
galaxies. Anyway, for our purposes, the behaviour at large radii is 
not that important, as the thermalization driven superwind solution 
is valid only close to the galax)lf| (see e.g. Strickland & Heckman 
2009). So, we will proceed to give the expression corresponding to 
the external gravitational potential. 

By taking the limit J? — > oo in equation l |13t , one finds that the 
total dynamical mass is given by Mt = Mdm(1 + Ai)^~°i. The 
expression of the associated gravitational potential for ^ a ^ 1 
is then similar to that given by equation (2) in Dehnen (1993): 



$h(ii) = 



Vcgjl + Aif- 
(2-ai)Ai 



R 



R + Ai 



2-ai 



(15) 



With this, we can establish new approximated thresholds for 
the open-box enrichment scenario and for accelerating superwind 
solutions. 



4 THRESHOLDS FOR OPEN-BOX ENRICHMENT AND 
ACCELERATING SUPERWIND SOLUTIONS 

When the effect of the external halo is considered, the asymptotic 
terminal speed is given by 



14 



5 - 2a 



\/2 



1/2 



(16) 



The flow enters into non- stationary regimes (inpouring or outpour- 
ing, see Paper I) when 



5 -2a 



y2 



(17) 



When the above inequality holds, the galaxy can eventually en- 
ter into an open-box enrichment scenario. Otherwise, we will 
have fully stationary solutions, unless radiative cooling or self- 
gravitation inhibit the stationary solution. 

Fully stationary superwinds have accelerating velocity profiles 
when 



-(2r,+ l):^$„(l) + 2r;fi±l 
V 5 - 2a 



(18) 



^ See also equation (3) in Dehnen (1993). 



^ This implies that the effect of the 'real' <I>h can be emulated there by 
giving adequate values to Ai and ai. 
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otherwise, they have decelerating velocity profiles. When the 
equality holds in the above relation, we have an almost constant 
external velocity profile with characteristic velocity 

= (2,,)-i/^[-2,^h(r,c)]^/^ = 



i) "rot 



(1 + 

(2-Qi)yli 



1 + 



2 — ai 



1/2 



(19) 



where v^^t is the rotation speed at r^, 



5 EFFECT ON THE HYDRODYNAMICS 

In our model, the dynamical mass (A/dm) contained within a bulge 
or galaxy nucleus experiencing an starburst episode is related to the 
concentration parameter and steepness of the external halo, and to 
the total dark matter and baryonic mass: 



Ml 



(20) 



(l + ^i; 

The dynamical mass contained in the external halo (r > rsc) 



Mh = Mt 



1 



(1 + ^1 



(21) 



Similarly, the dynamical mass contained up to an external 
characteristic radius (normalized to Tsc), Rt>, is 



Md = Mt 



(22) 



The radius _Rd can be associated to the 'disc' radius of the 
BM or well to the BM+DM virial radius. So, all the relevant galaxy 
parameters are correlated, in a similar fashion as in the work of 
Salucci et al. (2007). Nevertheless, we emphasize that the relation- 
ship between the parameters is alike but of course not the same, 
since here we constructed our theoretical model only following the 
results of the simulations of Abadi et al. (2010) and Navarro et al. 
(2010). 

We will proceed to discuss the effect of the extended haloes 
on the hydrodynamics. In order to do this, we consider the hydro- 
dynamical models presented in Table[T]and Table|2] The first table 
gives the inner parameters for three galaxies with different char- 
acteristics. The second table gives the properties of their external 
haloes. The groundwork for the discussion will be the premise that 
the spherical symmetric superwind solution is a zeroth-order ap- 
proximation to the aspherical case. We will consider again a refer- 
ence effective terminal speed of 2500 km s~^ for the case of null 
mass-loading, fully efficient thermalization, and total participation 
within the starburst volume, i.e. for e = /3 = = 1. For mod- 
els 2 and 3, the SFRs were obtained from formula (1) in Rupke, 
VeiUeux & Sanders (2005a) and formula (28) in Paper I, i.e. we 
considered SFRs that are consistent with the typical observed lumi- 
nosities for the object types, and that in parameter space, place the 
objects below the threshold for catastrophic cooling. We find that 
the predicted temperature profile is barely modified by the presence 
of the extended haloes. However, drastic changes are produced in 
the velocity profile. 

The first model is an extended version of model 5 in Paper 
I, and corresponds to a synthetic isolated dwarf elliptical galaxy 
that tries to emulate the characteristics of the most massive outlier 
of the mass-metallicity relationship detected by Peeples, Pogge & 
Stanek (2008, see also Paper I). We assumed that the galaxy formed 



by a contraction of ~ 40% of an initially unperturbed subhalo of 
DM and BM which had ^ 70% of its total mass located within 
r ~ 3rsc, so we used A\ — 0.5, Qi = 3/4 and Ru = 1. The latter 
is equivalent to saying that in this case there is no disc, i.e. we only 
have a galaxy nucleus. The internal dynamical mass distribution 
follows a plateau-like profile, which implies that some mechanism 
-perhaps the action of early powerful superwinds associated to a 
more extended and powerful starburst episode (see Governato et al. 
2010) or internal dynamical processes- has also transformed the 
initial mass configuration. In this model, starburst activity still per- 
sists near the galaxy centre, but with a high concentration. We as- 
sumed a low thermalization efficiency, which implies a small num- 
ber of massive stars and SNe within the characteristic concentration 
radius, A = 0.1. The justification for this is that the SFR is low, 
and that although small, the concentration radius is much larger 
than the typical radius of a massive star, i.e. the filling factor is 
low. Similarly, because of the small number of massive stars, just a 
small incorporation of mass is necessary to produce a heavily mass- 
loaded superwind. In this model, the presence of the extended halo 
suppresses the free superwind solution and the galaxy experiences 
an open-box enrichment [see equation ( |17H by keeping the met- 
als processed by the few massive stars still present near the galaxy 
centre. This will require however an already gas-poor galaxy at the 
moment at which the pollution occurred (Peeples, Pogge & Stanek, 
2008). As suggested above, the required low mass fraction could 
have been produced by the action of early superwinds associated 
to previous and more powerful starburst activity. This is consis- 
tent with the views of Peeples et al. (2008), which regarded their 
sample of outliers as transitional galaxies in their way to becoming 
typically isolated dE and dSph galaxies, but with a high metallic- 
ity. The suppression of the free superwind solution is practically 
insensitive to the value of ^5 ai 1, which indicates that the 
enrichment is produced by the physical conditions within A and 
the initial concentration of the unperturbed subhalo from which the 
galaxy formed. 

The second model considers the synthetic and very massive 
blue compact dwarf galaxy modeled in Paper I. Here, we add an 
extended disc to the model in order to 'transform' the galaxy into 
a luminous infrared on^fl (LIRG, Lm ~ 10^^ Lq). LIRGs and 
ULIRGs may be the end result of the merging of two moderate- 
size spiral galaxies and display traces of convergence to an ellip- 
tical morphology (Sanders & Mirabel 1996; Rupke et al. 2005a). 
We will model a LIRG that still exhibit evidence of heavily warped 
and thick discs, displaying a morphology perhaps similar to that 
of the central component of Arp 299 (Sargent & Scoville 1991; 
Heckman et al. 1999; Hibbard & Yun, 1999; Hu et al. 2004), but 
with just one nucleus. We assume that such features extend to up 
to 5 times the radius of the merger nucleus; thus, Rb = 5. The 
assumed mass and extension are consistent with CO emission ob- 
servations of (U)LIRGs (Lonsdale, Farrah & Smith 2006 and ref- 
erences therein). We further assume that the merging process has 
similarly transformed the steepnesses of the internal and external 
mass profiles of the interacting galaxies unperturbed haloes, such 
that a — Qfi = 3/4. We adopt the value Ai — 1 since it pro- 
duces interesting proportions. In such a case, ~ 66% of the BM 



N.B. As LIRGs and ULIRGs, BCDs may be the result of mergers, al- 
though generally they have lower masses, given that they mostly form from 
the merging of dwarf galaxies. Nevertheless, on a higher end, luminous blue 
compact galaxies can have dynamical masses of up to ~ 10" Mq (Gar- 
land et al. 2004, Pisano et al. 2010). 
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Figure 1. Superwind velocity profile for model 2 (solid line). The dashed- 
line represents the profile that would result if the external halo were ne- 
glected. 

and DM of both galaxies is contained within the warped discs char- 
acteristic radius and about ~ 30% of this fraction resides within 
the merger nucleus (that is ~ 20% of the total mass). As in the 
original model, starburst activity is present in the nucleus with a 
somewhat high concentration (A = 0.4), the thermalization effi- 
ciency is 0.5 and mass loading is important, /3 = 3 (see Heckman 
et al. 1999). In this model, the gravitational field of the external 
halo transforms the accelerating superwind solution associated to 
the original model into a bounded decelerating one (Fig. [TJ. This 
effect occurs because now we have a more massive galaxy. The 
produced deceleration will enhance the observable properties of the 
superwind because of a proportional density increment. However, 
an even larger total mass could inhibit the superwind solution. As 
mentioned in Paper 1, this is consistent with the superwind scaling 
properties found by Rupke, Veilleux & Sanders (2005b), whom re- 
ported and initial increment of the superwind observable properties 
with galaxy mass and a posterior flattening with the same. 

Rupke et al. (2005b) also reported a flattening of the super- 
wind observable properties at high SFR. In principle, the normal- 
ized free superwind solution (Section [2ll is insensitive to the SFR 
(provided that it could be considered constant during a relatively 
large time interval), as it just depends on the effective and asymp- 
totic terminal speeds. However, high SFRs will intensify the effect 
of radiative cooling, as more mass will be injected per unit time and 
volume, and thus, the stationary solution could also be radiatively 
inhibited. We have properly addressed this issue in section 2.1.2 of 
Paper I. 

As an extreme example of the effect of the nominal value of 
the galaxy mass, we model a massive and 'rare' radio galaxy with 
a very extended halo (see e.g. Genzel et al. 2003). We consider 
a galaxy with a dynamical mass of 1 x 10^^ M© within its nu- 
cleus of Tsc ~ 2 kpc. A mildly concentrated starburst (A — 0.5) is 
present in the nucleus, which has a cuspy dynamical mass distribu- 
tion (q = 1). We consider that the steepnesses of the inner region 
and the halo are the same and that the total mass of the galaxy is 
Mt = 4 X 10^^ Mq. This requires that Ai = 1. This implies that 
the half-mass radius is r ~ 2.5rsc and that ~ 80% of the total mass 
is contained within r ~ 8rsc- In this model a high deceleration of 
the superwind is produced, and the flow is unstable to small varia- 
tions of the effective terminal speed (thermalization efficiency), as 
shown in Fig [2] As a consequence, the flow could eventually en- 
ter into the outpouring or even the inpouring regime. On the other 
hand, if instead of a continuous steepness, we consider that the typ- 
ical cuspy halo profile (with slope a = 1) resulted from the con- 
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Figure 2. Velocity profiles for model 3. The solid lines corresponds to 
the parameters showed in Table [T] For this parameters, Voo ~ 1208 km 
s~^. The lower (upper) solid line (does not) consider(s) the presence of 
the external halo. Similarly, the dashed and dash-dotted lines con'espond to 
Voo = 1150 km s"'^ and Voo = 1125 km s~^, respectively. For the latter 
case, the stationary free superwind solution does not exist. 

traction of a smoother one, say with ai — 3/4 and Ai — 1, the 
free superwind solution would be inhibited and the galaxy could 
enrich itself with is produced metals in an open-box scenario. This 
would occur because in the second case, the total mass is slightly 
larger, Mt = 4.75 x lO'^^ Mq. The cumulative dynamical masses 
of the two assumed external profiles are very similar, their ratio 
varies from a value of 1 at rsc (they are identical as they must), up 
to a value ~ 0.86 at r = lOrsc; nevertheless, such a small varia- 
tion is enough to suppress the stationary superwind solution. This 
reflects the fact that at the limit of large galaxy masses, galaxies 
will retain most of their metals, as expected. 



6 SUMMARY AND CONCLUSIONS 

Here we have presented an analytic free superwind model that in- 
corporates the effect of extended DM and BM haloes. We find that 
the gravitational field of the extended haloes associated to massive 
galaxies can drastically alter the free superwind velocity profile and 
enhance its observable properties. We also find that massive haloes 
can also contribute to the inhibition of the superwind solution. 

In our model, the galaxy total mass (BM-fDM), the mass con- 
tained within a bulge or galaxy nucleus (defined by the character- 
istic radius rsc), the mass up to the disc characteristic radius, and 
the steepness and concentration of the external halo, are all corre- 
lated. Since the correlations are nonlinear, deviations from galaxy 
to galaxy are permitted, see Tables [T] and |2] Oppositely, there is 
no correlation between the above parameters and the concentration 
and steepness of the mass distribution for r < rsc. This is consis- 
tent with the results of the cosmological simulations carried out by 
Abadi et al. (2010) Navarro et al (2010), in the sense that haloes 
are not strictly universal. This should be expected, as we based our 
model in the 'structural contraction' property derived from their 
simulations. On the other hand, in their extensive work, Salucci 
et al. (2007) found that the previous parameters were correlated 
for spiral galaxies, and proposed universal rotation curves assum- 
ing a Burkert (1995) profile for the DM distribution. Our theoret- 
ical work diverges from theirs in that we considered additionally 
the mentioned inner concentration and steepness, which traces star- 
burst episodes. Such a consideration discards the possibility of uni- 
versal halo profiles and rotation curves, since in general they will 
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Table 1. Reference hydrodynamical models. Galaxy parameters for r < r. 



Model 


Type 


a 


A 


^SC 


Mdm 


SFR 


/3 


e 


c 


Vao 


Ve 


Regime (No halo) 










(kpc) 


(xlO* M0) 


Mq yr-i 








km s~^ 










(a) 


(b) 


(c) 


(d) 


(e) 


(f) 


(g) 


(h) 


(i) 


(j) 


(k) 


1 


dE 





0.1 


1 


100 


0.1 


4 


0.2 


1 


560 


0.2740 


borderline 


2 


(L)BCD/LIRG 


0.75 


0.4 


2 


500 


~ 40 


3 


0.5 


1 


1021 


0.2060 


accelerating 


3 


Radio 


1 


0.5 


2 


1000 


~ 200 


3 


0.7 


1 


1208 


0.2946 


accelerating 



Superwind hydrodynamical models. Table headers: (a) steepness parameter, (b) concentration pai'ameter, (c) radius, (d) dynamical mass, (e) star formation 
rate, (f) mass loading factor, (g) thermalization efficiency, (h) participation factor (i) effective terminal speed, (j) squared ratio of the escape velocity to the 
effective terminal speed, and (k) flow regime when the external halo is neglected. 



Table 2. Reference hydrodynamical models. External halo parameters. 



Model 


Type 


a\ 








Mt 


Regime 










(kpc) 


(XlO** Mq) 


(XlO** Mq) 








(a) 


(b) 


(c) 


(d) 




(f) 


1 


dE 


0.75 


0.5 


1 


A/dm 


~ 2.5AfDM 


open-box enrichment 


2 


LIRG 


0.75 


1 


5 


~ 3.16A/DM 


~ 4.76Mdm 


decelerating 


3 


Radio 


1 (3/4) 


1 


^hm = 2.5 


2A/dm 


4AfMD 


decelerating (open-box enrichment) 



External halo parameters for the models presented in Table[T] Table headers: (a) steepness parameter, (b) concentration parameter, (c) 'disc' radius (d) 'disc' 
mass, (e) total mass, and (f) Regime. In model 3, -Rhm corresponds to the half-mass radius. 



differ for r < r^c', however, the discrepancy will be reconciled at 
larger radii, and thus one could talk of an 'asymptotically universal' 
property, in the sense defined by Salucci et al. (2007). 

From the theoretical point of view, the importance of the an- 
alytic solution here presented resides in that it can be used to 
construct approximated superwind solutions when departures from 
spherical symmetry are important. In such a case, not just the pa- 
rameters of the galaxy but also its morphology will determine both 
the superwind hydrodynamics and the fate of the ejected gas; spe- 
cially, of metals. The connection of these features with the observed 
dispersion of the M-Z relationship will be discussed in a forthcom- 
ing work. 
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